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Abstract. The first successful multi-particle model for the host of the well-known quasi-stellar object (QSO) 3C 48 is reported. 
It shows that the morphology and the stellar velocity field of the 3C 48 host can be reproduced by the merger of two disk 
galaxies. The conditions of the interaction are similar to those used for interpreting the appearance of the "Antennae" (NGC 
4038/39) but seen from a different viewing angle. The model supports the controversial hypothesis that 3C 48A is the second 
nucleus of a merging galaxy, and it suggests a simple solution for the problem of the missing counter tidal tail. 
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1. Introduction 

3C 48 (Barkhouse & Hall 2001) was the first QSO to be dis- 
covered optically (Matthews et al. 1961; Matthews & Sandage 
1963) and the first QSO to be directly identified with a host 
galaxy (Boroson & Oke 1982). Its basic properties are listed 
in Table 1. 3C 48 has attracted much attention regarding the 
proposed evolutionary sequence of active nuclei (Sanders et al. 
1988). According to this scheme, interactions and mergers of 
galaxies trigger an evolution via ultra-luminous infrared galax- 
ies (ULIRGs) to QSOs. The observational evidence, however, 
is hampered by the fact that many transitionary objects show 
only dubious indications of past or recent mergers. Clarification 
requires detailed multi-particle modeling which helps with dis- 
entangling the complex spatial structure of merger remnants. 
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Table 1. Basic properties of 3C 48. 



RA 2 ooo 01h37m41.3s Ma et al. (1998) 

Dec 20 oo +33d09'35'.'l Ma et al. (1998) 

Redshift 0.367 Barkhouse & Hall (2001) 

Luminosity Distance 1,581 Mpc" 
Angular Size Distance 846 Mpc a 
Scaling Factor 1" « 4. 1 kpc 



" H = 75 km s 1 Mpc 1 and q = 0.5 will be adopted throughout 
the paper. 



3C 48 is an example of a transitionary object with pro- 
totypal properties in many respects: It has the typical far- 
infrared excess, originating from thermal radiation of dust 
which is heated by the quasar nucleus and by newly form- 
ing stars in the host galaxy (Neugebauer et al. 1985; Stockton 
& Ridgway 1991). Large amounts of molecular gas (Scoville 
et al. 1993; Wink et al. 1997) indicate the possibility of a young 
stellar population in the host. Finally, long-slit spectroscopy 
gives evidence for an ongoing starburst in the host which cur- 
rently seems to be close to its maximum activity (Canalizo & 
Stockton 2000). But the merger scenario for 3C 48 is still un- 
clear: Indeed, the host has a significant tail-like extension to the 
northwest whose tidal origin is rather compelling with regard 
to the kinematics and ages of its stars (Canalizo & Stockton 
2000). However, the nature of the apparent second nucleus 
3C 48A about 1" northeast of the QSO (Stockton & Ridgway 
1991) and the location of the expected counter tidal tail remain 
an unsolved problem. 3C 48 A could as well be due to the radio 
jet (Wilkinson et al. 1991) interacting with the dense interstellar 
medium (Chatzichristou et al. 1999). A feature at the southeast 
of the 3C 48 host, previously interpreted as a counter tidal tail 
(Boyce et al. 1999), has turned out to be a background galaxy 
(Canalizo & Stockton 2000). Instead, a counter tail extending 
from the southeast to the southwest is suspected (Canalizo & 
Stockton 2000) but not yet identified. Canalizo & Stockton 
(2000) suggest that such a location of the two tidal tails might 
be explicable by a certain projection of the merger scenario 
used to simulate the "Antennae". 

This paper reports the first successful multi-particle model 
for the 3C 48 host. Suggesting simple solutions for the 3C 48A 
problem and the counter tail problem, the model largely re- 
solves doubts about the merger hypothesis for 3C 48. 
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2. Methods 

2.1. Numerical simulations 

The stellar-dynamical 3-dimensional N-body simulations are 
performed with TREESPH (Hernquist & Katz 1989), a tree 
code used in its non-collisional mode. Stable initial particle dis- 
tributions for the model galaxies are set up with BUILDGAL 
(see Hernquist 1993). The spatial orientation of the two disk 
galaxies is parameterized by their inclinations i with respect 
to the orbital plane and their pericentric arguments cj as intro- 
duced by Toomre & Toomre (1972). Both galaxies have ex- 
tended mass distributions so that the orbit of their encounter is 
not Keplerian but decaying. Two descriptions can be used to 
characterize these orbits: The first one is a pseudo-Keplerian 
description using the parameters of eccentricity e, pericentric 
distance r per i, and angle to pericentre £2 per i of the correspond- 
ing Keplerian orbit for which the total mass of each galaxy is 
associated with a point mass at the respective centre of mass. 
The second one is a direct description of the decaying orbit us- 
ing the true apocentric and pericentric distances r apo and r per ; 
of the first passage to define an eccentricity in its generalized 
formulation e = (r apo - r per i)/(r apo + r per i). For convenience, 
the system of units, which remains intrinsically scale-free, is 
scaled to the system suitable for 3C 48. 

The results of the simulations are analyzed as 2- 
dimensional projections. In order to mock the pixel array data 
of imaging observations, the particles are sorted into a 5 12x5 12 
grid. The virtual pixel values are computed by adding up all 
particles located in a grid cell along the line-of-sight. Without 
any special weighting of a nuclear component, the mock im- 
ages are comparable to QSO-subtracted images of the 3C 48 
host. Spectra for each grid cell are generated by sorting the par- 
ticles into velocity channels according to their respective line- 
of-sight velocities. Thus, an average stellar line-of-sight veloc- 
ity is assigned to each virtual pixel. The resulting data arrays 
are spatially smoothed by Gaussian convolution and converted 
into FITS format to facilitate the further data processing with 
standard astronomical software. 



2.2. Observational data on 3C 48 

The data presented by Canalizo & Stockton (2000) are used 
for comparing the simulations with observations. They provide 
information about the optical surface-brightness of the QSO- 
subtracted 3C 48 host (Fig. 1 therein) and about the stellar 
kinematics along the four slits A, B, C, G (Fig. 1 and Table 2 
therein). In reference to these data, the basic proportions of the 
main body of the 3C 48 host are classified by dimensionless 
length ratios (left panel of Fig. 1 and left column of Table 3). 
Such a comparison is independent of the length scaling, in 
contrast to the comparison of line-of-sight velocities which re- 
quires a positioning of the four slits on the mock image. Having 
determined the final physical length unit of the simulations, the 
curvature of the northwestern tidal tail is compared by using an 
angle-versus-distance plot (right panel of Fig. 1 and Fig. 5). 




Fig. 1. Sketch of the parameters used for characterizing the dimen- 
sions of the 3C 48 host. The lengths (LI, L2) are measured along 
the longest extension of the host body (PA « 4° for the observations, 
PA ss 19° for the simulations), the widths (Wl, W2) are measured per- 
pendicular to this. The curvature of the tail is traced along maximum 
intensity (right panel) by angle-distance-vectors (a,, rf ; )- 

Table 2. Initial parameters of the two identical galaxies in the sys- 
tem of units suitable for 3C 48. Each galaxy consists of a spherical 
non-rotating bulge, a rotating exponential disk, and an isothermal halo 
(Hernquist 1993). 



Parameter 


Bulge 


Disk 


Halo 


number of particles 


8,000 


8,000 


8,000 


softening length [kpc] 


0.21 


0.28 


1.4 


mass [10 10 M o ] 


1.86 


5.60 


32.48 


scale length [kpc] 


0.88 


3.50 


35.0 


maximum radius [kpc] 


7.0 


52.5 


105.0 


scale height [kpc] 




0.7 




maximum height [kpc] 




7.0 





3. The 3C 48 look-alike 

Different mass ratios of the initial galaxies, different snapshots 
during the merger process, and different projection angles of 
the merger remnants were probed in a still limited parameter 
study. 

The nearest 3C 48 look-alike is found for the merger of two 
identical galaxies whose physical and numerical properties are 
given in Table 2. With these parameters the galaxies are similar 
to spirals of type Sb. The experimental setup is the same as 
used for simulations of the "Antennae" - i.e. both galaxies are 
symmetrically oriented with i\ = z'2 = 60° and <x>\ — u>2 — 
-30° (see Toomre & Toomre 1972; Barnes 1988). The model 
galaxies are initialized near the apocentre of the corresponding 
elliptical Keplerian orbit which is defined by the eccentricity 
of e = 0.5, the pericentric distance of r per ; = 20 kpc, and the 
period of 1 .2 Gyr. The time step in the simulations is fixed to 
1 Myr which guarantees that energy is conserved to better than 
1% during the merger. The true decaying orbit is characterized 
by the generalized eccentricity of e = 0.8 and the pericentric 
distance of r per i = 7 kpc. The outer regions of the two galaxy 
bulges begin to merge after ~ 272.5 Myr, just before pericentric 
passage. The merging is not a single process but the centres 
of the two bulges are repeatedly flung apart before they settle 
down in one common density peak. 
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Fig. 2. Two different projections of the same simulation snapshot after 461.1 Myr. The left panel shows the projection for which the merger 
remnant looks like the "Antennae", the right panel shows the projection for which the remnant looks like the 3C 48 host. The coordinate planes 
indicate the respective tilt of the orbital plane (x-y). See text for a detailed description. 
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Fig. 4. Comparison of the observed and simulated line-of-sight velocities vlos along the four slits A, B, C, G. Offsets are measured in the 
direction indicated by the arrows in Fig. 3. Dots represent the average velocities of the 3C 48 look-alike. The area enclosed by v LO s + Av L os 
and v L0S - Avlos (dashed histograms) corresponds to the confidence region of the stellar line-of-sight velocities measured for the 3C 48 host. 
These data are taken from Fig. 1 and Table 2 in Canalizo & Stockton (2000). 
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Table 3. Measured proportions of the 3C 48 host, as taken from the 
contour plot in Fig. 1 of Canalizo & Stockton (2000), compared to the 
proportions of the nearest look-alike. 



Ratios 


3C 48 host 


Look — alike 


L1/L2 


0.76 


0.70 


Ll/Wl 


1.25 


1.44 


L1/W2 


1.67 


1.86 


Ll/D 


0.52 


0.55 


L2/W1 


1.65 


2.06 


L2/W2 


2.20 


2.64 


L2/D 


0.69 


0.79 


WI/W2 


1.33 


1.29 


Wl/D 


0.42 


0.38 


W2/D 


0.31 


0.30 
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Fig. 3. Contour plot of the surface brightness of the 3C 48 look-alike. 
Length units are fixed in arcsec by positioning the four slits A, B, 
C, G used by Canalizo & Stockton (2000). The small inset shows a 
magnified view on the still separated density peaks of the merging 
bulges. 



The nearest 3C 48 look-alike emerges after 46 1 . 1 Myr. Two 
projections of this merger remnant are shown in Fig. 2. In the 
left panel ("Antennae" look-alike), the view is perpendicular to 
the orbital plane x-y. In the right panel, (3C 48 look-alike) the 
orbital plane is tilted southwards, westwards, and counterclock- 
wise by 120°, 160°, and 116°, respectively. The proportions of 
the 3C 48 look-alike are listed in the right column of Table 3. 
The positions of the four slits A, B, C, G on the look-alike 
host and the resulting physical coordinate system are shown in 
Fig. 3. In Fig. 4, the scaled velocities along the slits are com- 
pared to the confidence region of stellar line-of-sight velocities 
given for 3C 48 by Canalizo & Stockton (2000). The angle- 
versus-distance comparison for the curvature of the northwest- 
ern tidal tail is shown in Fig. 5. 
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Fig. 5. Comparison of the curvature of the northwestern tidal tail of 
the 3C 48 host (solid line) and its look-alike (dashed line). 



4. Discussion 

4.1. The nature of 3C 48 A 

Optical and near-infrared images of the 3C 48 host show two 
luminosity peaks at the positions of 3C 48 and 3C 48A, the 
latter being located about 1" northeast of 3C 48 (e.g. Stockton 
& Ridgway 1991; Canalizo & Stockton 2000). With softening 
lengths of ~ 0.25 kpc (0706) for the bulge and disk particles 
(see Table 2), the spatial resolution of the simulations is high 
enough to identify corresponding density peaks in the 3C 48 
model. As shown in the small inset in Fig. 3, the centres of 
the bulge components of the two merging galaxies are still sep- 
arated at the stage of the 3C 48 look-alike. Their distance of 
about 076 (2.5 kpc) and their relative positions on a southwest- 
ern to northeastern axis are similar to the observed configura- 
tion of 3C 48 and 3C 48A. Thus, a scenario with 3C 48 and 
3C 48A being the two centres of merging galaxies is possible. 
However, the exact configuration of the density peaks in the 
simulations is very sensitive to the projection angle and to the 
time at which the snapshot is taken. About 20 Myr later, the two 
peaks have already merged into one. Furthermore, the purely 
stellar-dynamical model does not address the question of a pos- 
sible nuclear activity at the positions of 3C 48 and 3C 48A. A 
detailed discussion of nuclear activity depends on whether or 
not a black hole exists at the mentioned positions and on the 
respective fueling rates. 

4.2. Location of the counter tidal tail 

Since each of the two merging disk galaxies forms a tidal tail, 
the missing second tidal tail has always been a caveat of the 
merger hypothesis for 3C 48. Here, the simulations suggest a 
simple solution: At the projection angle of the 3C 48 look-alike, 
the second tidal tail is mainly located in front of the body of the 
host and, therefore, severely foreshortened. It extends from the 
southwest towards the northeast, roughly along slit B, so that 
measurements along this slit trace a mixture of line-of-sight 
velocities from the body and from the tail. This could explain 
why the observed and the simulated line-of-sight velocity sig- 
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natures along slit B (Fig. 4) are dominated by scattering around 
a mean velocity close to zero. Slits A, C, and G, in contrast, 
are characterized by large absolute line-of-sight velocities (up 
to ~ 200 to 300 km s~') and strong variations along the slits. 
A counter tidal tail extending from the southwest towards the 
northeast in front of the main body of the 3C 48 host is a com- 
pletely new alternative. Regarding the information about stellar 
kinematics, this location seems to be more likely than the two 
suggested tails in the southeast and from the southeast towards 
the southwest (Boyce et al. 1999; Canalizo & Stockton 2000) 
which have failed identification so far. 
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4.3. The evolutionary history of 3C 48 

Conclusions about the orbital parameters for 3C 48 and the 
original parameters of the merging galaxies can only be ten- 
tative. The orbital period of the best fit model amounts to 
about 20% of the age of the universe at the redshift of 3C 
48 (~ 5.4 Gyr). A merger scenario with such an orbital pe- 
riod is plausible, assuming an initially highly eccentric orbit 
of the merging galaxies which is transformed into a bound or- 
bit by dynamical friction of their dark matter halos (e.g. Jones 
& Stein 1989). It has been found that the morphology and the 
kinematics of tidal tails are very sensitive to the rotation curve 
of the interacting model galaxies (Dubinski et al. 1996; Mihos 
et al. 1998; Dubinski et al. 1999). Thus, instead of two iden- 
tical galaxies, an alternative model for 3C 48 could start from 
two galaxies with different rotation curves so that only one of 
them forms an extended tidal tail. However, even in its general- 
ity the multi -particle model presented in this paper gives rather 
compelling evidence that the formation of 3C 48 is linked to 
a merger process. Therewith, 3C 48 ranks among these transi- 
tional objects which support the evolutionary scenario (Sanders 
et al. 1988) in its original merger-driven definition. 
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